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In an accelerated exclusion process (AEP), each particle can "hop" to its adjacent site if empty 
as well as "kick" the frontmost particle when joining a cluster of size I < £ max - With various choices 
of the interaction range, ^ ma x, we find that the steady state of AEP can be found in a homogeneous 
phase with augmented currents (AC) or a segregated phase with holes moving at unit velocity (UV). 
Here we present a detailed study on the emergence of the novel phases, from two perspectives: the 
AEP and a mass transport process (MTP). In the latter picture, the system in the UV phase is 
composed of a condensate in coexistence with a fluid, while the transition from AC to UV can 
be regarded as condensation. Using Monte Carlo simulations, exact results for special cases, and 
analytic methods in a mean field approach (within the MTP), we focus on steady state currents and 
cluster sizes. Excellent agreement between data and theory is found, providing an insightful picture 
for understanding this model system. 

PACS numbers: 64.60.De 64.75.Gh 05.60.-k 89.75.Fb 



I. INTRODUCTION 

Unraveling rich behaviors emerging from simple ingre- 
dients in systems driven far from equilibrium is a con- 
tinuous pursuit in theoretical physics. Non-trivial flux 
of physical quantities, or current, captures the macro- 
scopic feature of the complex system and is often closely 
governed by the intrinsic dynamics. There are numerous 
examples where the steady state current depends sensi- 
tively on the constituents in the system through both 
long-range and short-range interactions such as queueing 
in traffic and pedestrians [1], transport of biomolecules 
[2H5] and minerals j6]. 

One of the venerable models, the totally asymmet- 
ric simple exclusion process (TASEP) not only provides 
many interesting mathematically exact results of non- 
equilibrium statistical mechanics [TrfTT] , it also brings in- 
sights to potentially important applications in, for ex- 
ample, protein synthesis [12f[15| and regulating vehicu- 
lar traffic [16] . Closely related is the zero-range process 
(ZRP) in which particles hop between sites with rates 
determined by the origin site occupancy. The mapping 
between ZRP and an ordinary TASEP proved illuminat- 
ing in our later discussions. ZRP also finds its versatility 
in studying granular materials as well as phase separation 
in one-dimensional systems |17j . 

Inspired by assisted hopping in transcription [21 H] , we 
introduced a new variant of TASEP recently [IS], the 
"accelerated exclusion process" (AEP), to characterize 
interactions beyond nearest neighbors among particles. 
Violating detailed balance, AEP is non-equilibrium in 
nature and involves a few surprising phenomena. Let us 
briefly recapitulate the dynamic rules of TASEP before 
turning to the novel features of AEP: In each update at- 
tempt of TASEP, a particle in a discrete one-dimensional 



(ID) lattice of L sites is chosen at random to hop into 
its neighboring site (provided that is vacant) with rate 
7 (typically chosen as unity). Both periodic and open 
boundary conditions (particles enter and exit with rates 
a and /3 respectively) have been studied extensively [3- 
lllj . The non-equilibrium steady states of an ordinary 
TASEP are well-understood. For example, the current- 
density relationship is given by Jtasep = pO- — p), m the 
thermodynamic limit. 

In the AEP with periodic boundary condition, N par- 
ticles are placed in a ring of L sites and follow the same 
rules as in TASEP. In addition, when a particle hops to 
a cluster of particles of size £ < ^ max , it simultaneously 
"kicks" the frontmost particle of that cluster one site for- 
ward. There is no avalanche, as the "kicked" particle does 
not trigger another kick. When ^ max = 0, AEP reduces 
to an ordinary TASEP. A schematic of AEP is shown in 
Fig. 0a). 

In this article, we focus our attention on the non- 
equilibrium steady states of this AEP and explore the 
interplay amongst the overall density p = N/L, £ max and 
current J(p;Anax)- Due to "kicking," a number of novel 
features arise in the AEP. Given that the contribution 
to the current can be either 1 (a hop) or 2 (a hop and 
a kick) [TB], the system is naturally expected to display 
an augmented current (AC). In this state, the system is 
homogeneous and we may anticipate that Jac lies be- 
tween Jtasep and 2 x Jtasep- It is intriguing that, at 
high densities (p > 1/2), Jac can exceed 2 x Jtasep- 
Meanwhile, for low densities, the kicking action results 
in "facilitated" or "cooperative motion" [TM2"T] . where 
the average velocity of the particles, v, can be increased 
by adding particles to the system. Clearly absent in 
the ordinary TASEP, this phenomenon can be charac- 
terized by d p v > or a positive curvature in J: d 2 p J > 
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FIG. 1: (a) AEP with £ max = 4. A particle can hop (hollow 
arrows) to its unoccupied right neighboring site and kick (grey 
arrow) a second particle when joining a cluster of size I < 
Lax. (b) Mapping AEP in (a) to MTP. The ball in stack 3, 
if chosen, lands immediately in stack 5 (hop and kick), while 
the one in stack 1 lands in stack 2 (hop). 



[19H2T] . Furthermore, for moderate values of £ max and 
at high densities, the system exhibits an inhomogeneous 
state in which the particles "condense" into a macro- 
scopic, "solid" cluster, in coexistence with a "fluid" of 
density ~ 1/2. Surprisingly, the (average) current is just 
1 — p = H/L, where H is the number of holes in the 
system. Thus, the fluid can be regarded as a loosely 
bound set of holes, moving together with unit velocity 
(UV). These two different states of the system will be 
referred to as the AC and the UV phases, respectively. 
As the overall density p is increased with fixed £ max , we 
observe a discontinuous jump in J(p;£ ma ,x), from a non- 
trivial Jac(p) t° the simple Jyy = 1 — p. Our goal is to 
understand these remarkable phenomena. 

Although AEP is originally cast in the language of an 
exclusion process, both the intuitive picture and the anal- 
ysis for predicting the aforementioned features turn out 
to be much easier when viewed in an equivalent represen- 
tation, the mass transport process (MTP). In the next 
Section, we present a detailed description of the AEP- 
MTP mapping, which is a simple generalization of the 
TASEP-ZRP mapping. In this setting, the exact master 
equation can be easily written. Following a brief sum- 
mary of simulation results in Section [Till we provide the- 
oretical considerations for the properties of AC and UV 
phases in Sections IV and[Vj respectively. In Section VI 



configuration of the system, C, can be characterized by 
the set of site occupancies, {n,}, i = 1, L, with n being 
or 1. Alternatively and more conveniently, AEP can be 
formulated as a mass transport process (MTP) [52] in 
which we regard the particles in front of each hole as 
"balls" stored in a "stack." Each hole in AEP becomes 
a stack in MTP, labeled by a = 1, ...H. The £ a balls in 
stack a correspond to the cluster of particles between the 
a th and (a + l) th hole. Each stack may be occupied by 
any number of (indistinguishable) balls: £ a = 0, 1, N. 
An equivalent specification of C can thus be the set {£ a } 
instead. 

The mapping from AEP in Fig. fija) to MTP is shown 
in Fig.[TJb). The dynamic rules of AEP thus become: In 
each update attempt, a random stack a is chosen. If it 
is empty (l a = 0), another attempt is made. Otherwise, 
one of the balls in the chosen stack hops to the next 
stack and if £ a +i G [l,^max], the ball takes a second hop 
immediately, landing in stack (a + 2). For instance, the 
"hop and kick" scenario in Fig. [lja) becomes a hop from 
stack 3 to 5 in Fig. [IJb). When £ max = 0, it returns to 
an ordinary TASEP. 

Clearly, N = T, a £ a is conserved. These rules are sum- 
marized in the master equation for P(C;t) in Eq. ([TJ). 
One of the advantages of MTP representation is that {£ a } 
provides directly the cluster size distribution in the AEP, 
a measure which we can use to characterize the AC and 
UV phases quantitatively. 

The only non-trivial aspect of this mapping is J, the 
overall particle current. In each update attempt, a ran- 
dom site is chosen for AEP, while in MTP is a stack. 
Thus, the time scale differs by a factor of H/L. In partic- 
ular, if we compute the average number of particle move- 
ments in MTP, it must be multiplied by H/L = 1 — p 
when compared to J in AEP. 

We proceed to formulating the dynamic rules as a mas- 
ter equation. To facilitate this task, we define the char- 
acteristic functions 



1 if£ e [1,4 



otherwise. 



we venture a phase diagram for this system. We conclude 
and provide an outlook for further quests in Section [Vll| 



II. ACCELERATED EXCLUSION AS A MASS 
TRANSPORT PROCESS 

Regarded as particles traversing a ID ring, AEP allows 
particle to move only when it is adjacent to a hole. The 



Note that £ max is an implicit parameter in these func- 
tions. The master equation governing the evolution of 
P(C;t), namely the probability to find the system in 
configuration C = {£ a } in t attempts after some initial 
configuration, for H > 2 is: 
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X {8 (£' a+1 ,£ a+ l) 6 (4+2,4+2 + l) *(4+l) + S (4+1,^ + 1 + 1) 5 (4+2,4+2) *(4+l)} 



Here, 5 is the Kronecker delta. For H = 1, there is just 
one stack in MTP, leaving AEP trivial. For H = 2, the 
problem is easily solvable and hints at the significant role 
of tax- The next case (H = 3) is the first non-trivial one 
and the exact solution in the L — ¥ oo,£ max — > oo limit 
provides valuable insight into the UV phase. We show 
the details of H = 2 and 3 in Appendices [A] and [B| 

The configuration space {£ a } consists of the lattice 
points in an (H — l)-dimensional hyper-tetrahedron. The 
easiest way to visualize this is its standard embedding in 
-H-dimension, i.e., a (linear) space joining the following 
H points: (N, 0, 0), (0, 0, N). For H = 3 and 4, 
the configuration space is just an equilateral triangle and 
the standard regular tetrahedron, respectively. 

In general, this dynamics does not obey detailed bal- 
ance, so that finding an explicit stationary distribution, 
P* (C), is not simple [23, 24J. To have some understand- 
ing of the system behavior, we will exploit the approxima- 
tion schemes presented below, before which let us com- 
ment briefly on some general properties of our system. 

Two extreme cases are noteworthy. One is £ m&K = 0, 
which is simply the ordinary TASEP. Dropping the \ 
term, Eq.Q reduces to 

1 H 

^tasep (C; t + 1) = Jj J2 Ptasep ( C; *) x 

c a =i (2) 

x<*(4,4- i)5(4+i,4+i + i) n Hfatp) 

Though this dynamics also violates detailed balance, it 
does satisfy the "pairwise balance" condition |7j, so that 
■Ptasep (^) ^ ■*-• The opposite extreme is £ max — > oo, or 
simply £ max > L or N. x(£) reduces to 5 (£) so that {...} 
in Eq. becomes: 

{<5(4 +1 ,4+i) 5(4+2,4+2 + i)£(4+i,o) + 

+5 (4+i,4+i + l)5(4 +2, 4+2 5(4+1,0)} 

where 5=1 — 5. This important limit exemplifies the 
AC phase and will be examined more closely below. 

Returning to the general case, one can compute the av- 
erage of any quantity Q in the stationary state assuming 
P* is known: 

where the sum is taken over {ri;} or {£ a }, whichever is 
more convenient. In particular, the average particle cur- 



rent in steady state is given by: 

J (p; 4ax) = (rii (1 - n 2 ) (1 + n 3 )) 

- /ni(l- n 2 ) f[ n 3+m \ (4) 

\ m=0 / 

= (l-p) (5(4,0) [1 + X(4)]) 

Note that we have invoked translational invariance of 
P* in writing these expressions. If the reader is concerned 
that this invariance may be (spontaneously) broken, as 
in the UV phase, then J = £ {n} (...) /L = (...) /H 

may be used instead. Also, the second line is the result of 
having one (two) hop(s) for a particle-hole-hole (particle- 
hole-particle) triplet. In MTP, whenever stack 1 is cho- 
sen, a ball will hop, with an additional hop if the target 
stack has the requisite number of balls. 

III. SIMULATION RESULTS 

We exploit the random sequential updating scheme and 
simulate AEP with a typical lattice of L = 1000. Initial- 
izing the system with N particles on successive lattice 
sites, we make N attempts to update in every Monte 
Carlo step (MCS) so that each particle has on average 
one chance to hop. The contribution to the current from 
both hops and kicks are instantly accounted for. The 
typical length of our simulations is r = 5 x 10 5 MCS. We 
start data collection after 2 x 10 5 MCS to ensure that 
the system has reached steady state and all measure- 
ments are averaged over 3 x 10 5 MCS thereafter unless 
otherwise specified. 

First let us examine the AC phase. We note that not 
only is J always larger than Jtasep except at p = and 
1, it even surpasses 2 Jtasep when p > 1/2, as shown in 
Fig.[2j The simulation (red circles) is for £ max — 1000 and 
thus all hops lead to kicks. The system is always in the 
AC phase regardless of the density, which may lead one to 
believe that Jac should be close to 2 Jtasep = 2/9(1 — /?), 
since a typical particle will hop two sites instead of just 
one. This rough estimate is reasonably good for high 
densities, as Fig. [2] (black dot-dash line) indicates. Using 
the most naive mean field approximation, we can start 
from the exact expression Eq. Q and let £ max — > oo. 
This leads us to p (l — p 2 ). Surprisingly, this provides a 
much poorer overall picture (blue dotted line in Figs. [2] 
and [3| . Furthermore, there is some inflection at low den- 
sity where the curvature of Jac(p) is positive, indicating 
an "accelerated region." This is illustrated in Fig. [3] 
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FIG. 2: (Color online) Comparison amongst Jac from 
Monte Carlo simulation, Jtasep, 2Jtasep, a naive mean 
field approximation and an improved mean field approxima- 
tion through MTP (MTP-MF). The result from MTP-MF in 
Eq. (Ill (blue solid line) provides remarkable agreement with 
the simulation. 
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FIG. 3: (Color online) Curvature of J(p) from simulation 
(circles), a naive mean field approximation (blue dotted line) 
and MTP-MF (blue solid line). 



Utilizing the MTP picture, we provide a mean field the- 
ory (MTP-MF) in Section [TV] to describe the AC phase, 
which yields exceptional agreement with simulation data, 
shown in solid blue line in Figs. [2] and [3] This approach 
leads to a remarkably good description of all phenomena 
presented here. Moreover, it provides a viable explana- 
tion for why the estimate p (l — p 2 ) fails. The details are 
presented in Section |TV| 

The emergence of the UV phase depends on the ap- 
propriate choice of ^ max . In all of our simulations where 
L = 1000, there exists a UV phase for 10 < £ max < 950. 
We show a few typical £ max 's in Fig. [4] When the system 
enters UV, J becomes independent of € max and is simply 



(1 — p), indicating the holes are moving at average speed 
1. It is therefore both more effective and illuminating to 
switch to the reference frame of moving holes: Starting 
from a full system with one giant cluster of particles, i.e. 
H = 0, as more and more holes are injected to the sys- 
tem, we find the system settles into a cluster of particles 
in one part and a half-filled region in the other part of 
the system. Exploiting both AEP and MTP frameworks, 
we explicate the emergence of such "phase separation" in 
Section [V] 




FIG. 4: (Color online) Transition from AC to UV for ^ max = 
20 (red circles), 100 (green squares) and 500 (yellow triangles) 



IV. AUGMENTED CURRENT PHASE: 
FIELD DESCRIPTION 



MEAN 



In this section, we focus on the stationary state of the 
AC branch, in which the system is homogeneous and the 
density profile is just p. Now, if we let £ max > L, then 
the system be in AC for all p's. In the MTP represen- 
tation, the rule is especially straightforward: Every hop 
into an occupied stack makes a further hop. Thus, the 
probability that a stack is occupied 

f=(S(tx,0)) 

will play a central role. Further, we find it useful to study 
the more detailed distribution 

P* = 

i.e., the probability that a stack contains precisely i par- 
ticles. Of course, / = X^ >0 P* (-0 implies 



1 - / = P* (0) 



A. Steady state occupations 



(5) 



Although it is not possible to derive an exact master 
equation for P(l;t) = £ 8 {£, l{) P (C; t) from Eq. 0, 
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we will exploit a mean field approach to find P* (£). To 
maintain a steady state, a given stack must gain and lose 
a ball with equal probability in any attempt. In other 
words, we can find P* by balancing the average rates of 
gain and loss. 

Clearly, the probability that an occupied stack (£ > 
1) loses a ball is P* (£) (with the trivial factor of 1/H 
suppressed). Meanwhile, the stack can gain one ball (i.e., 
I— 1 — > t) if both stacks upstream are occupied. Injecting 
the mean field approximation, we estimate this condition 
by f 2 . If the chosen stack is empty, there is an additional 
way it can gain, from its immediate upstream neighbor 
if occupied. Thus, the balance equations read: 

(/ 2 + /) P* (0) = P* (1) 

fP*(i) = P*(£+1);£>0. 

Strictly, we should account for the upper limit £ < N. If 
we neglect such finite size size effects, the normalization 
constraint is £^ P* (£) = 1, leading us to an explicit 
expression for the stationary distribution: 

P*(£>0) = (l + /)/ 2f - 1 P*(0) (6) 
P*(0) = 1-f (7) 

Note that Eq. ^ is entirely consistent with Eq. ^5b. 

Next, let us relate / to the density p. Computing (£) 
is straightforward: // (l — f 2 ). But this average occu- 
pation must be N/H = p/ (1 — p), leading to 



(S) 



Instead of a cumbersome algebraic expression for f(p), 
let us define r\ by 



sinh 7/ 



H 1 -p 



2N 2p 

and recognize (l — f 2 ) // = / _1 — /, so that 



(9) 



(10) 



or explicitly, / = exp < — sinh 



1 



2p 



1311 . Not sur- 



prisingly, / is a monotonically increasing function of p. 

In the next subsection, we exploit these results to find 
the current-density relationship and explore some conse- 
quences. 



B. Currents and velocities 

In MTP, an occupied stack (/ > 0) contributes one (or 
two) hop depending on whether the next stack is empty 
(or occupied), the probability of which is associated with 
(1 — /) (or /). Within our approximate scheme, the to- 
tal contribution is / [(1 - /) + 2/] = / (1 + /). Together 



with the H/ L factor to scale from MTP to AEP, we con- 
clude that the mean field approximation for the particle 
current is 



tMTP-MF 

j ac 



(p) = (l-p)/(l + /) 



(11) 



This simple approximation provides an excellent predic- 
tion, with no fit parameters. The solid blues line in 
Figs. [2] and [3] show its remarkable agreement with our 
simulation data. 

Recall p [1 — p 2 ), the "naive mean field approxima- 
tion" for J (p), which is far from simulation data (see 
Fig. [2] except near p = 0,1). What is the difference be- 
tween the two approaches - that one performs so much 
better? Comparing Eq. (Ill with p (l — p 2 ), we see 
the difference to be simply J[l + f) instead of p (1 + p). 
Since / is an estimate of the probability of successive 
holes not being nearest neighbors, the failure of p here im- 
plies that it underestimates much of the clustering of the 
particles,. The remedy provided by MTP can be cast as a 
simple and intuitive picture, for which we coin the phrase 
"abhorrence of empty stacks." Since an empty stack can 
be filled by hops from two (upstream) stacks, it is far less 
likely to remain empty, compared to the situation in the 
ordinary TASEP. It is easy to check that f — pis positive 
for p e (0,1) and peaks at ~ 0.12 (around p ~ 0.57). 
In other words, / is more successful at accounting for 
the scarcity of empty stacks, consecutive holes in AEP 
language (at moderate densities). 

Now let us consider "velocity." There are two notions 
of velocity associated with J (/?). One is the average par- 
ticle velocity, v = J/p. From Eq. ( 11 ), we obtain a simple 
expression: 



.MTP-MF 



J/p=(l-f 2 )(l + f) 



(12) 



Since particles never hop backwards, this v is necessarily 
positive, approaching 1 (or 0) in the limit p — > (or 1). 

The other notion of velocity is d p J. Similar to the 
group velocity for waves, d p J is sensitive to collective be- 
havior such as the motion of fluctuations or disturbances. 
Thus, it is negative at high densities, corresponding to 
holes moving "backwards" (the simplest case being the 
single-hole system). An unusual and notable feature of 
AEP is "cooperative motion" [IT)l421) . In an ordinary 
TASEP, adding a particle to the system always reduces 
both v and d p J . In AEP, by contrast, there is a regime in 
which particles "cooperate" and move faster when more 
are present: Namely both d p v and d 2 J can be positive. 
Both simulation data and d f v MTp - MF = (1 + /) (1 - 3/) 
display such a regime. Similarly, this behavior is also 
present in d 2 J, studied in [TM2"T] and easily discerned as 
curvature shown in Fig. [3j Given the agreement in Fig. [2] 
it is not surprising that the MTP-MF is also successful 
at predicting the phenomenon of cooperative motion. 

Additionally, let us comment on two other features as- 
sociated with the current J^T P_MF (p). It is clear from 
Eq. ^ that, at p = 0.5, we have N = H and so 
/(! + /) = 1. Inserting this result into Eq. (Ill, we 
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conclude that J^JP-mf ( Q ^ _ Q 5 Thig ig larger th&n 
simulation data by about 1%, a difference which may 
be attributed to either statistical errors in the data or 
MTP-MF's failure to capture certain correlations. Since 
our system does not obey explicit particle-hole symmetry, 
it is unclear if J = 0.5 is merely a curious coincidence or 
an indication of some hidden symmetry. Of course, this 
is the "maximal" value of 2</tasep, but again, this co- 
incidence may also be accidental. Finally, we note that 
both the data and J^J P_MF exceeds 2Jtasep for much 
of p > 1/2, both peaking around p = 0.54. Apart from 
the mathematical analysis, we have not developed a good 
intuitive picture for explaining these observations. 



Cluster size distribution 



Exploiting the MTP-MF further, we can compare the 
result in Eq. ^ with measurements of the cluster size 
distribution (CSD). In Fig. [5} we show simulation data 
for N = 600 and 900, both with £ max = 1000. It is clear 
that the distributions are consistent with exponentials. 
Further, if we construct ratios of the successive values, 
we find the average over a range (where the scatter of 
the data is small) to be 0.556 and 0.895, respectively. For 
comparison, using Eqs. (|9 10 1 to compute / 2 , we find the 
values being 0.519 and 0.895, respectively. Another point 
for comparison is P* (0), for which the simulation data 
provide 0.291 and 0.0582 versus 0.279 and 0.0540 from 
Eqs. ([7j [9j [To]), respectively. Such agreement leads us 
to conclude that the MTP-MF approach indeed captures 
the essence of the AC phase. 
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FIG. 5: (Color online) Cluster size distributions in AC phase. 
C — L — 1000 with N = 600 (magenta diamond) and 900 
(black triangle), respectively. Theoretical predictions are not 
shown explicitly for clarity. The agreement between them and 
data are comparable to those in Figs. [2] and [3] 



V. UNIT VELOCITY PHASE: PERSPECTIVES 
FROM BOTH AEP AND MTP 

In this section, we focus on the UV branch, in which 
the system is inhomogeneous and exhibits an approxi- 
mately half-filled region in coexistence with a fully occu- 
pied domain. The presence of this phase depends cru- 
cially on having a moderate £ max . The most remarkable 
feature of this phase is that the average particle current is 
exactly (1 — p), shown in Fig. [4] In other words, the aver- 
age speed of the holes (or the "disturbances" ) is precisely 
1, independent of particles being added to, or removed 
from, the system (until a phase boundary is reached). 
In both AEP and MTP representations, there exist sim- 
ple descriptions which provide an intuitive and appealing 
picture for why such an unusual state can persist. The 
following subsections are devoted to each of these per- 
spectives. Before presenting the details, note that such 
a state spontaneously breaks the translational symmetry 
underpinning the dynamics. Interestingly, this broken 
symmetry is manifested in slightly different forms in the 
two representations. We will comment on this difference 
in each of the subsections below. 



A. A "hole train" in AEP 

In the language of the original exclusion process, this 
phase is best understood if we focus on how the holes 
move. When we choose the particle at site k — 1 and 
find that it can hop to a hole at site k, we can regard 
this process as choosing the hole and exchanging it with 
the partner particle. Next, we should ask if there is a 
(particle) cluster of length I occupying the sites from fc+1 
on. When I £ [l,i max ], we also move the hole at k+£ + 1 
to k + 1. In other words, when a hole moves, it will pull 
the next hole ( "downstream" ) , provided the gap between 
them lies in [l,4nax]- We should remind the reader that 
the second hole does not pull a third one. 

Since this phase is present only in the high density 
regime, it is natural to first consider systems with the 
lowest values of H. These provide us the necessary pic- 
ture to understand the existence of a UV phase. A system 
with a single hole evolves trivially: The hole "swims up- 
stream" with unit average velocity, since every attempt 
to move it is successful. The H — 2 case is somewhat 
more interesting. Deferring the details to Appendix [XJ 
we state the principal results here. It is straightforward 
to enumerate, for any L and ^ max , all possible station- 
ary states, through which the important role played by 
4nax is revealed. Most significantly, for moderate £ max /L, 
the two holes form a tightly bound pair. To be precise, 
the gap between them can be either or 1, with equal 
probability. Let us emphasize that this is an absorbing 
state. Starting with any initial separation, the two holes 
will eventually drift together, form the bound state, and 
never become unbound thereafter. Since the "leading" 
hole always moves when it is chosen, the pair moves with 



7 



unit average velocity. 

The first non-trivial case is H — 3. For a finite system 
with arbitrary ^ max , an exact solution is not yet avail- 
able. Nevertheless, we can gain some insight into the UV 
behavior through an exact result in the limit L — > oo 
followed by t max — > oo (details in Appendix |B|) . Here, 
the leading pair remains tightly bound, as the third hole 
cannot affect the leading hole. Denoting the gap between 
the second and the third hole by m, we know that the 
third hole can trail behind the second by m > 0, so that 
the complete description of the system lies in the follow- 
ing probabilities. Let po,i (w) be the probability that the 
gap between the first hole pair is or 1, with the third 
hole trailing by m sites. We illustrate examples of p (5) 
and pi (4) in Fig. [6j As shown in Appendix |Bj we find 
(apart from m = 0, 1) 



Po,i {m) oc C r 



(13) 



where £ = 2 — \[2 = 0.586. In other words, the third 
hole is "loosely" bound, with an exponential tail of a 
characteristic length /i = — 1/ (ln£) = 1.87. The intuitive 
picture is clear: The first pair advances with UV, with 
the third hole being "pulled along," at a typical distance 
fi behind the second. It is natural to label such a triplet 
a "hole train" with the first pair being the "engine" [TS] . 

(a)po(5) 



n nnnnn n 



(b)p.(4) 



FIG. 6: Sketch of a "hole-train" with H = 3: (a) Two leading 
holes are consecutive with the third trailing by m = 5. The 
probability of this scenario is denoted as po(5). (b) Two lead- 
ing holes have one particle in between, with the third trailing 
by m = 4, denoted as pi(5). 

Continuing this line of thought, we see that, in the 
limit posed above, the first pair will always form an "en- 
gine" which advances with UV, unaffected by how many 
holes trail behind them. While we do not have exact so- 
lutions for the general case, we can argue why the holes 
should be bound and the image of a train is quite reason- 
able. In particular, consider the last hole - the "caboose" 
in the parlance of freight trains - and the penultimate 
hole, which we will refer to as "X." If X is chosen and 
moves, it will pull the caboose along (provided the gap 
in-between is non-zero). The gap remains the same. On 
the other hand, if the caboose moves, the gap decreases. 
Only when X is pulled along by the hole upstream from it 
does the gap increase. Although these considerations ap- 
pear to imply that the gap between penultimate and final 
holes performs a typical random walk, we point out that 
the entire train length does not increase even when this 



gap increases. Thus, we argue that this mechanism can 
"hold the train together." Accepting this scenario, we 
see that adding or removing holes to the system merely 
changes the total length of the train. Meanwhile, since 
the engine advances with UV, the whole train also moves 
as such, leading to an i7-independent velocity. In the 
following section, we present a more tenable and quan- 
titative argument for the existence of this hole train, as 
well as an estimate of its average density (~ 0.5) in the 
language of MTP. 

Before proceeding, let us comment on two other as- 
pects of the AEP perspective. First, although a typical 
snapshot of our system clearly violates translational in- 
variance, this symmetry is restored quite quickly (O (L) 
MCS) since the hole-train moves at UV. To be precise, 
if we measure the occupation at a specific site, it will 
settle at p within such times. If, on the other hand, we 
tag a particular hole and measure the occupation in one 
of its nearest neighbor sites, then the results will expose 
the inhomogeneity inherent in the system. Restoration of 
the symmetry (in finite systems) would take much longer 
than O (L) MCS, a subject well beyond the scope of this 
work. Second, in a finite periodic lattice, the average dis- 
tance between the caboose and the engine, A, is finite. 
For systems with A ^> f miK , the caboose does not af- 
fect the (lead hole in the) engine. Thus, the integrity 
of the engine remains intact and moves the entire train 
with UV. As holes are added or removed, A becomes 
smaller or larger, respectively. When enough holes are 
added, or if £ max is raised, then A can approach £ max , 
the caboose can destroy the engine, and the train can 
become unbound. In a nutshell, this is the mechanism 
for the transition from UV to AC, as the overall density 
is lowered (e.g., Fig. [4]). We will return to this picture in 
Section IVH 



B. 



Condensation in MTP 



As for the AC phase, the MTP representation provides 
us with a more quantitative picture. The solid parti- 
cle cluster plays the role of the condensate in non-trivial 
ZRP's. With a finite (and small enough) £ max , it is pos- 
sible for one stack to contain more balls than £ max . Such 
a stack can gain balls in two possible ways, much like 
how an empty stack can be filled above. Meanwhile all 
stacks can lose a ball in just one way. Thus, the number 
in this stack will grow, until a steady state is reached. 
It is natural to refer to this behavior as "condensation" 
and this stack (or the balls in this stack) as the "con- 
densate." To make contact with the previous subsection, 
note that there are A balls in the condensate stack. Need- 
less to say, each of the remaining H — 1 stacks is likely 
to hold very few balls. In the language of ZRP, the state 
of the other stacks is referred to as "a fluid." Here, we 
recognize them as the hole train. Further, this picture 
allows us to appreciate better why the hole train remains 
bound. First, the train length is monotonically related 



to the fluid density, becoming longer/shorter when the 
condensate loses/gains balls. Second, the only way to 
redress the imbalance (two gains vs. one loss) for the 
condensate is when the fluid density remains relatively 
low, i.e., a set of stacks with few balls in each. This 
scenario corresponds to a bound train. 

Turning to a more quantitative description of the 
steady state, we denote the occupation probability within 
the train by /train- A good estimate for it comes from the 
balance of the gain/loss contributions from the conden- 
sate, namely 



/train + /train 1 



(14) 



The predicted value, /train = (a/5 — l) /2, is not very 
illuminating. Instead, by comparing with Eq. ([8]), we 
find a more insightful relation: 



H 



train 



= N t , 



(15) 



namely, p t rain = 0.5, implying a train length is 27J train . 
This result also indicates that the typical distance from 
one hole to the next in the train is around 2, a picture 
entirely consistent with the result \x ~ 2 in the H = 3 
case. 

Meanwhile, since the train consists of H — 1 stacks, we 
arrive at -Htrain — H. Further, we have At ra in + A = N = 
L — H so that the size of the condensate is given by: 



A = L - 2H = L (2p - 1) 



(16) 



All these predictions are borne out relatively well in 
simulations. For example, in Fig. [7j we show the CSD's 

= 100. Clearly, the con- 



fer N = 600 and 900 with L 



800 



densate sizes are seen to fluctuate around 200 and 
respectively, as predicted by Eq. (16). The properties of 



the fluid/hole-train, as revealed by the small clusters not 
shown in Fig. [7J are essentially the same in both cases. 
For small clusters the distribution indeed decays expo- 
nentially, with ratios of the successive values being 0.406 
and 0.412, respectively. These are somewhat higher than 
/ t 2 rain = 0.382. Given that our theory is based on a mean 
field approximation, we speculate that the difference are 
due to non-trivial correlations, the study of which is be- 
yond the scope of this work. 

So far, the analysis is focused on the average behav- 
ior of the fluid and the condensate. Since our approach 
considers the single-stack occupation, we can apply it to 
the condensate and exploit a self-consistent way to pre- 
dict P* on (A), the probability for the condensate to have 
A balls in the steady state. Note that, unlike A, A is a 
variable here. In such a configuration, the fluid has only 
A— A balls, which allows us to estimate / (the occupation 
probabilit y of a stack in the fluid) as a function of A. Us- 
ing Eqs. (9 10), we find / (A) = cxp j-sinh -1 2( ^~ X A) j 

while / (A) is just /train- We can now use / to estimate 
the rate at which the fluid supplies a ball to the conden- 
sate, namely, / 2 + /. Denoting this rate by 



we find an expression similar to Eq. ([6| 

5 (A)P c * on (A) = P c * on (A + l), (18) 

since the condensate loses at unit rate. It is straightfor- 
ward to check that g (A) is a monotonically decreasing 
function and is unity at A — A. Thus, P* on (A + 1) = 
P* on (A) = P are the peak values of the distribution and 
can be conveniently used to start the recursive evalua- 
tion of two sequences: P* on (A + 1 + k) and P* on (A - k) 
with k > 0. Furthermore, even though g appears to de- 
pend on both control parameters (H, N), it actually is a 
function of a single (shifted and scaled) variable 



A-A A-N+H 



H - 1 



H- 1 



(19) 



For completeness, we provide the explicit expression: 

g(A;H,N) = 1 - ^ exp j- sum' 1 ^ -L j (20) 

Now, we can express In P c * on as a sum over In g : 

lnP* n (A + l + fc) = ln£-£> fl (0 (21) 



lnP c * on (A-fc) = In ^ + 5^ In 5 (0 (22) 



where the sums run over £ being integer multiples of 
1/ (H — 1), up to ±fe. Although we cannot evaluate this 
sum, we can extract its properties for large H and k of 
O (1) (i.e., small £'s). To leading order, the resultant 
lnP c * on is a function of k 2 /H. In other words, the con- 
densate size distribution (at this level of approximation) 
is universal in the following sense: Although its explicit 
dependence is P* on (A; N, H), it can be cast in scaling 
form, P c * on cx $ (x), where $ is a universal function of 
the scaled variable 



A-N + H 



H 



(23) 



(17) 



Such behavior is similar to Gaussian distributions, which 
are universal apart from a displacement and a rescaling. 
A detailed study of $ is in progress and will be reported 
elsewhere. Here, let us present the numerical results from 
Eqs. plp2| for the cases above. The agreement with the 
two data sets are again remarkably good (Fig{7j). We note 
the slight discrepancies in the A = 200 case, and believe 
that they are the consequences of the fluid section being 
longer (H being 400 instead of 100). Surely, for systems 
with larger fluid components, the fluctuations therein will 
be more serious. Obviously a careful study and analysis 
of such fluctuations and correlations will be necessary if 
the goal is to go beyond mean field theory. 

To summarize, we see that a quantitatively coherent 
picture of the hole train emerges when viewed in the MTP 
representation. Here, the hole-train corresponds to the 
fluid, while the solid cluster in the rest of the lattice cor- 
responds to the condensate. Translational symmetry in 
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(b) 
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FIG. 7: (Color online) Cluster size distributions in UV phase. 
Cm = 100, L = 1000. Circles are simulation data and the 
solid line is the result from MTP-MF. (a) N = 600. (b) 
N = 900. 



the MTP is spontaneously broken, as the condensate re- 
sides in one of the H stacks. However, unlike in the AEP 
picture, this condensate does not move with UV from 
stack to stack. Symmetry restoration must proceed by 
evaporation and re-condensation. Such a process is ex- 
pected to take considerably longer than the O (L) MCS 
in the AEP representation, while its detailed nature is 
being investigated [25] . The results here allow us to take 
the thermodynamics limit: L, N, H, £ max — > oo with fi- 
nite 



P = N/L- X = £ max /L. 



(24) 



Provided A < 2/9—1, it is possible for a hole train to form, 
occupying a finite fraction (2 — 2p) of the ring. The gap 
between the engine and the caboose, corresponding to 
the size of the condensate, fills the remaining fraction: 
2/9—1. This result clearly implies that the UV phase 
cannot exist for p < 1/2. 



VI. TRANSITIONS BETWEEN AC AND UV 

In this section, let us consider the transition between 
the two phases and map out a phase diagram in the 
P ~ ^max plane. First, note that the thermodynamic 
limit cannot be studied rigorously, especially since the 
exact steady state distribution, P* (C), is not known. If 
this limit does not exist, then the standard term "phase" 
should be used with some caution. Second, the standard 
approach to phase transitions involves taking this limit 
with the stationary state, i.e., the limit t — ¥ oo is taken 
first, while simulations are based on running finite sys- 
tems (L < oo) for finite times (r < oo). Therefore, we 
can only make some estimates and offer some rough ar- 
guments here. Obviously, more convincing conclusions 
can be drawn from a thorough finite-size scaling anal- 
ysis (See e.g. a task beyond the scope of this pa- 
per. Finally, we should comment on the order parameter, 



i.e., how we characterize the phases. By studying mainly 
the current J, we have implicitly chosen (the operators 
in) Eq. Q here. Yet, as discussed in the previous two 
sections, the phases may be better characterized by the 
presence/absence of a macroscopic cluster. Thus, an- 
other possibility is to study the distribution of the size of 
the largest cluster, which we denote by Q (s). Deep in the 
AC p hase, Q should be similar to P* (£) for large I. From 



Eqs. ( |6|10| ), we therefore expect Q (s) — > e 2r,s . On the 



other hand, deep in UV, this cluster is the caboose-engine 
gap or the condensate, so that Q (s) is just P* on (A). In- 
deed, most of our theoretical arguments for the phase 
transition presented below will be based on the proper- 
ties of Q. 

Since we are dealing with non-equilibrium steady state, 
there is no widely accepted notion of a "free energy" even 
if we managed to find an explicit P* (C). Thus, we can- 
not follow the standard route, defining a first order tran- 
sition through a jump in its derivative. Alternatively, we 
can define such a point dynamically, given that our sys- 
tem is formulated as a stochastic process. A reasonable 
choice is, for example, that set of control parameters with 
which the system settles for equally long periods in each 
of the phases while switching between them occasionally 
("tunneling" ). If computer power/time is unlimited, we 
can measure s (t) for arbitrarily long periods and com- 
pile a histogram for Q (s). If our analytic power is strong 
enough, we can access the exact steady state Q. For a 
range of parameters, Q should be sharply bimodal, al- 
lowing us to locate special points where the modes are 
equally probable. However, as both approaches are quite 
limited at present, we can offer only rough estimates and 
reasonable arguments for a "phase diagram" here. Three 
regimes are expected to be present: pure AC, pure UV, 
and "mixed" (AC+UV). The best way to characterize 
these regimes is through lifetimes. Specifically, deep in 
the pure regimes, the system settles relatively quickly 
into one phase, regardless of initial conditions. By con- 
trast, in the mixed regime, once it settles into AC or UV 
(typically through judicious choice of initial conditions), 
that state can persist for extraordinarily long times. In 
particular, it is possible for these lifetimes to scale expo- 
nentially with the system size, L. In that case, such a 
regime rightly deserves the label "bistable" [32] . 

We emphasize that there are three independent control 
parameters in the simple AEP. They can be L, p, and A, 
or N, H, and £ max , for example. Also, a variety of "ther- 
modynamic limits" can be taken, depending on the order 
that different quantities are sent to infinity. In addition, 
for simulations, the initial condition and length of runs 
will be important to consider. Enumerating all possibil- 
ities is exceedingly difficult, if not impossible. Here, we 
will focus mainly on the parameters used in our simula- 
tions (L = 1000 and a wide range of p and ^ max ) and 
provide some arguments from our analysis for other sit- 
uations. 

It is clear that if we use a totally inhomogeneous sys- 
tem (all particles clustered together) as an initial con- 
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dition, then we are likely to find the pure UV regime, 
as well as to explore the boundary between the mixed 
and pure AC regimes. Indeed, all the simulation data 
presented above are collected under these conditions. As 
indicated in the previous section, in AEP, UV is associ- 
ated with a "hole-train" of typical length 2H, which will 
be destroyed if the caboose wanders within ^ max of the 
engine. Thus, UV is unstable if € max > L — 2H . This 
provides an estimate for the critical density associated 
with the boundary with the pure AC phase |18j : 



L(2p cA -l) 



(25) 



To be explicit, a system with p < p C A = will set- 
tle into AC only. This boundary is shown as the dashed 
(blue) line in Fig. [8] (with L = 1000). Although the agree- 
ment with data (red circles) is reasonably acceptable, this 
estimate can be improved by incorporating some fluctu- 
ations. In MTP, the caboose-engine gap appears as the 
condensate size, A, which fluctuates around A. Thus, A 
can reach £ max with a small probability, P* on (^ max ), even 
if A may not be near £ max . Now, we may argue that, in 
a run of r MCS, rare events with probability 1/r can oc- 
cur. Exploiting this connection and approximating the 
scaling form for P* on (A) by a Gaussian, we see that the 

nax ~ O (VFhi7 



UV can become unstable if A 



Inserting r ~ 10 used in our simulations, we find that 



Eq. ( 25 ) is slightly modified. As an illustration, we plot 



L (2p cA - 1) - y/L{l - p cA ) lnr. 



(26) 



Shown as the solid (blue) line in Fig. [8j it is arguably an 
improvement. If this result is upheld in a more rigorous 
analysis, we may conclude that if the thermodynamic 
limit is taken first, there is a non-trivial region in the p-X 
plane associated with AC+UV bistability, while Eq. ( 25 ) 
marks its border with the pure AC regime. 

Next, we explore the stability of the AC phase. In 
simulations, this phase will be more favored by distribut- 
ing particles uniformly on the lattice initially. In this 
manner, we expect to find another boundary, beyond 
which the system never settles in AC. Deferring a sys- 
tematic investigation, we simply provide a few examples 
in Fig. [8] (red diamonds). Theoretically, our approach is 
similar to the one above: What is the probability for the 
largest cluster to reach ^ max particles? Within our ap- 
proxi mate sc heme, this is given by P* (^ max ), and using 
Eqs. ( 6|9|10 ), we arrive at exp (— 2?7^ max ). Applying the 
connection to runs of length r, we obtain 



lnr 



(27) 



where r\ c u is related to the critical density associated with 
the boundary UV regime in Eq. ([9]): 



sinh rjcu 



PcU 



2p 



•cU 



(28) 



The resultant is also plotted in Fig. [8] (dot-dash blue line). 
While the discrepancies between this estimate and data 



is larger than those above, we may conclude that this 
approach is a viable first step. In particular, we believe 
that a major difference between these two cases lies in 
the following. For a system in UV to tunnel to AC, the 
condensate stack needs to wander all the way down to 
< mal . We can formulate this problem of finding the life- 
time as a first passage time of a single walker arriving 
at a particular destination. By contrast, in the reverse 
process, tunneling from AC to UV requires only one of 
the stacks to wander up to £ max , corresponding to finding 
the first time that any one of the H walkers arrives at the 
destination. Clearly, the latter problem is more complex, 
especially since the walkers are not entirely independent. 
To improve on Eq. (27) will be a worthy next step. 
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FIG. 8: (Color online) Phase diagram in p-^ max plane. Sym- 
bols are from s imulations. Phase boundaries are results from 
Eqs. ( 25|26|28 1 (dash, solid and dot-dash), respectively. 



To summarize, we presented a plausible phase diagram 
associated with the discontinuous transitions observed in 
simulations, consisting of three regimes. In two of these 
regimes, the system appears to evolve to a unique steady 
state: AC or UV. In between, our simulations show that 
the system can settle into either state, depending on, 
e.g., initial conditions. We conjecture that our system 
supports the phenomenon of bistability, namely, the time 
scales for switching between these states grow exponen- 
tially with they system size. In other words, we expect 
the behavior here to resemble that in equilibrium systems 
with long range interactions (See., e.g., [27] and more re- 
cently [351 US])- To prove or disprove this conjecture will 
likely be accomplished through careful observations of 
hysteresis along with a finite size scaling analysis. 



VII. SUMMARY AND OUTLOOK 

In this article, we investigated an accelerated exclusion 
process (AEP) on a ring where particles hop when the 
neighboring site is empty, as well as kick another one for- 
ward when joining a cluster of particles of size t < f max . 
Through Monte Carlo simulations, we discovered that, 
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with various choices of density p and interaction range 
< max , the system may be found in an augmented current 
phase or a unit-velocity phase. The behavior this AEP 
exhihibits, both dynamic and static, are much richer than 
the standard TASEP. Focusing on the steady state, we 
expand the findings reported in Ref. [18| and seek a com- 
prehensive theoretical framework for understanding these 
novel features. The apparent inadequacy of a naive mean 
field approach prompted us to seek alternative routes. 
Treating AEP as a mass transport process (MTP) of 
balls contained in stacks and may jump either one stack 
( "hop" only ) or two ( "hop and kick" ), we provide a 
more intuitive picture of both phases and the transition 
between them. In this representation, a mean field ap- 
proximation scheme is formulated to compute several key 
quantities. With no fit parameters, the predictions agree 
remarkably well with results from simulations. 

For the AC phase, we found an expression for the par- 
ticle current, J^JP-mf = (i _ p )f(\ + /) ) w ith / being 
the probability of an occupied stack (equivalently, the 
frequency of isolated holes in AEP), given explicitly by 
Eqs. ( 9|10 ). This result enabled us to quantitatively es- 
timate both J and the "acceleration" in the facilitated 
region. Additionally, it provided an intuitive picture for 
the scarcity of hole pairs, which leads to more "kicks" 
and augmented currents. 

Once the system is in the UV phase, J is simply (1 — p) 
regardless of ^ max , indicating the holes in the system 
are traveling at unit velocity. This intriguing result can 
be appreciated from the AEP and MTP representations 
with different insights. In the language of AEP, the sys- 
tem in UV is composed of a "hole-train" of length 2H, 
led by an "engine" composed of a tightly bound hole-pair. 
The complement of the hole train is a cluster of particles 
of size A = L - 2H. In the language of the MTP, the 
hole train and cluster is, respectively, the fluid and the 
condensate. We also computed the leading term in P* on , 
which enabled us to understand the average sizes of the 
condensate as well as its fluctuations. Preliminary stud- 
ies of scaling behavior and a universal distribution are 
encouraging and further investigations are in progress. 

Deferring charting a precise phase diagram the AC and 
UV phases to our further quests, we reported the essen- 
tials in the formation of condensates, thus infer the phase 
boundary using p and £ max as order parameters. Starting 
as a "solid" (all balls in one stack) or a "liquid" (balls 
distributed through all stacks) leads the system to favor 
UV or AC. Various factors, including the initial condi- 
tions, affect where the system eventually settles and how 
long it remains, hence we conjecture the two different 
phase boundaries presented in Fig. [8] with support from 
our simulations. 

Although this study provided valuable insights into the 
AEP, there are many avenues to improve on both sim- 
ulations and theory, in order to advance a better un- 
derstanding of its behavior. Examples mentioned above 
include explorations of the dependence on L and r, fi- 
nite size scaling analysis, and a careful study of clusters' 



evolution and the size distributions. On the theoreti- 
cal front, we should account for some correlations in the 
system, for example, by considering the joint distribution 
P (^i 7 ^2) in the MTP. A more refined phase diagram than 
our Fig. [8] would be most desirable. In particular, we may 
expect that the discontinuous jumps give way to a contin- 
uous, second-order like, phase transition. Subsequently, 
all the standard issues associated with such a transition 
can be explored, from critical exponents and universal- 
ity classes to scaling and renormalization group analyses. 
Beyond static properties, we envisage many interesting 
dynamic questions. In addition to investigations already 
in progress [35] , it would be instructive to study time se- 
ries and power spectra of various quantities, since they 
can expose the details of correlations in time. For exam- 
ple, we may study microscopic currents associated with 
entry and exit times from the condensate. While the 
latter is expected to be simply Poisson distributed, the 
former may be more complex, as it is connected to the 
fluctuations of the entire fluid. In particular, the correla- 
tions of entry/exit times should provide information on 
propagation of fluctuations (through the fluid) . 

Beyond the system studied here, there are natural gen- 
eralizations, such as having two or more particles being 
activated and ^-dependent kicking probabilities. Another 
natural generalization is the AEP with open boundary 
conditions, with a variety of injection/extraction possi- 
bilities. Mapping out the equivalent of the open TASEP 
phase diagram fully will be an arduous, but rewarding 
task. We may introduce inhomogeneous hopping rates 
modeling blockages or adsorption/desorption along the 
entire chain. Further afield, we may wish to consider 
systems with many species, or many lanes ("quasi- ID" ), 
as well as in higher dimensions. Yet another important 
task for the future is to see to what extent the features 
of AEP are present in more realistic models of systems 
in nature that display assisted hopping. Finally, we hope 
that AEP will be a new window for understanding not 
only exclusion processes, but also non-equilibrium statis- 
tical mechanics in general. 
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Appendix A: Exact solution for H — 2 

For H = 2, the system is sufficiently trivial that we 
can simply enumerate all possibilities, denoted as a pair 
of integers in parenthesis. In the MTP representation, we 
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need to consider only the number of balls in one stack, 
£G [0, N]. The other stack contains N — £. Furthermore, 
there are at most three intervals in [0, N], in which the 
rules are different. 

If tax > N, then every ball moves two steps (return- 
ing to the original stack), so that every "interior" (i.e., 
£ G [1,N — 1]) configuration stays the same. Meanwhile, 
each of the two "boundary" configurations decays as 2 _t , 
since choosing the filled stack will lead to a stationary 
one. Though it appears to be stationary in the MTP rep- 
resentation, the AEP current is always 2 no matter which 
particle-hole pair is exchanged. It is natural, therefore, 
for us to give such a state the label AC. Since every ini- 
tial condition corresponds to such a state, the system is 
always in AC. 

If N > £ max > N/2, then we can have a non-maximal 
"interior region" (N — £ max < £ < £ maK ) of stationary 
configurations, as in the previous paragraph. If the ini- 
tial condition is in this region, the system is again AC. 
However, if the system starts in one of the two "boundary 
regions," then it will evolve as follows. Such a configu- 
ration consists of one stack with £ > ^ max balls and the 
other with N — I < £ max . If the latter is occupied, then a 
ball leaving the first stack will make two hops and return 
to the original stack. Thus, this stack either gains a ball 
or remains the same, so that £ tends to drift upwards. In 
other words, the system performs a biased random walk 
towards the boundary until it reaches the configuration: 
(N, 0). From there, it can only reach (N — 1, 1). From 
this point, the system jumps between these two configu- 
rations, i.e., a stationary state we recognize as the tightly 
bound pair (the "engine") discussed in Section [V] A. It is 
clear (and straightforward to prove) that such a system 
should be labeled by UV. Thus, a system with initial 
conditions in these "boundary regions" simply evolves 
towards a UV stationary state. 

Finally, if N/2 > £ max , the configurations in the "inte- 
rior region" (£ metX < £ < N — £ max ) consist of both stacks 
having more than £ max balls. Now, balls just move from 
one stack to the other, so that the system performs an 
unbiased random walk in £. When it reaches one of the 
"boundary regions," it converts to performing a biased 
random walk as above. Thus, such system will always 
end in a UV stationary state. 

To summarize, the H = 2 case, though seemingly triv- 
ial, provides the essentials of the AC and UV "phases." 
The simplest "phase diagram" emerges: A domain in N- 
£ max plane with pure AC, one with pure UV, as well as 
a third where both AC and UV can be the end state 
(depending on initial conditions). 



Appendix B: Solution for H = 3 in a special limit 

Clearly there are more possibilities for the H = 3 case 
and they would be more complicated. Enumerating them 
and providing exact solutions in each scenario remain 
to be completed. Here, let us consider a special limit, 



L — > oo followed by £ max —¥ oo. The former limit implies 
that the caboose cannot affect the engine, which remains 
a tightly bound pair. The latter ensures that the engine 
can affect the caboose, which can lag behind by an ar- 
bitrary number (to > 0) of sites. In this scenario, we 
only need to consider po (to) and pi (to) , the probability 
that the gap between the first pair is and 1, respec- 
tively, with the caboose trailing by to (see Fig{6]) . In the 
stationary state, these satisfy: 



2p (to) = po (to + 1) + pi (to) (1 - 5 m0 ) 
+ 8 m opo (0) + 6 

HI 1P1 (o) 
3pi (to) = po (m) + pi (to + 1) 

+ pi (to - 1) (1 - <5 m0 ) + S m0 pi (0) 



(Bl) 



A simplification occurs when we sum the two sets (bal- 
ancing the currents between the po's and the pi's): 



Po M = Y V\ M = \ 



(B2) 



m>0 



m>0 



the last "=" being the result of normalization. To obtain 
the individual p's, we consider the generating functions: 



G. (z) = Z ™P' ( m ) 

m>0 



(B3) 



and verify that they satisfy 



1-2* z \( G Q \ = f (1- z) [zpi (0) + po (0)] \ 
1 l-z)\G x ) \ p„(0) + (l + z)pi(Q) J' 



where the second line expresses the balance of total fluxes 
between to and to + 1. Thus, we have, e.g., 



Go 



(1 - zf [po (0) + z Pl (0)] - z [pq (0) + (1 + z) gi (0)] 

2{z-z){z- l/2z) 

(B4) 

where 2 = 1— 1/V2. Since Go cannot be singular at 
2 < 1, the numerator must vanish there and leads to a 
relation between p (0) and p\ (0): 

(1 - zf [po (0) + z Pl (0))] - z [po (0) + (1 + z) Pl (0)] = 

(B5) 

A second relation between them comes from Eq. (|B2|) . 



1/2 = Go (1) = Po (0) + 2p x (0) 
and allows us to arrive at: 



Po (0) 
Pi (0) 



11 - 6^ 
14 

3^2-2 ^ 
14 ^ 



a 0.1796 
0.1602. 



(B6) 

(B7) 
(B8) 



Instead of writing explicit expressions for all the p's, let 
us exploit a shortcut to the asymptotic behavior, namely, 



subtract ( B5 ) from the numerator in ( B4 1 and cancel the 
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(z — z) in the denominator. The result is that Go (z) 
must be of the form 



Go (z) 



A + Bz + Cz 2 
1 - 2zz 



(B9) 



where A, B, C are constants that can be explicitly com- 
puted. From here we find that, for all m > 2, 

Po(m)=C m {A + B/( + C/t 2 } (BIO) 

where £ = 2z = 2 — v2- A similar expression can be 
derived for pi(m). Thus we see that the third hole is 
bound with an exponential tail, of characteristic length 
-l/ln£« 1.8697. 

The exact average of m, the distance between the en- 
gine and the caboose, can also be calculated via 

= d[G {z) + Gl (z)] u=i 



dz 
9%/2 



14 



1.4806, 



(Bll) 



a result slightly smaller than the characteristic length of 
the exponential tail, — l/ln£. 

As a cross check, we compute the average current ex- 
plicitly, via 

Po (0)+2 Po (l)+2p (2) + ...+3 Pl (0)+5 Pl (l)+5 Pl (2) + ... 
which is 

2 Po M + 5 Pi M - Po (0) - 2pi (0) = 3 



m— 



m=0 



by Eqs. (B2 B6). Since iJ = 3, this result shows the UV 



property explicitly. 

Needless to say, if we reverse the order of limits 
(< max — > oo followed by L oo), then the system will be 
only in the AC phase, even though the exact P* is yet to 
be obtained explicitly. 
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